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Abstract 



An alternative to the standard cold dark matter model has been recently 
proposed in which a significant fraction of the energy density of the universe 
is due to a dynamical scalar field (Q) whose effective equation-of-state differs 
from that of matter, radiation or cosmological constant (A). In this paper, 
we determine how the Q-component modifies the primordial inflation gravita- 
tional wave (tensor metric) contribution to the cosmic microwave background 
anisotropy and, thereby, one of the key tests of inflation. 

PACS number(s): 98.80.-k,95.35.+d,98.70.Vc,98.80.Cq,4.30.-w 



1 



Typeset using REVT^X 



I. INTRODUCTION 



A key prediction of inflationary cosmology is a nearly scale-invariant spectrum of initial 
fluctuations composed of energy density (scalar metric)0B§EI and gravitational wave (tensor 
metric)!' !0 perturbations. The two components leave imprints on the cosmic microwave 
background (CMB) anisotropy that are potentially detectable and distinguishable. The 
spectrum of each component is characterized by a spectral index (n) which determines how 
the perturbation amplitude varies with wavelength. Each spectrum is also predicted to be 
nearly scale- invariant; using the standard notational convention, this corresponds to n$ ~ 1 
for the scalar metric fluctuations and rix ~ for the tensor metric fluctuations. The spectral 
indices and the ratio of tensor to scalar amplitudes on any given scale are determined by 
the equation-of-state during inflation or, equivalently, the inflaton potential. A prediction of 
inflation is a simple relation between ns, nr and the ratio of tensor to scalar amplitudes that 
applies to nearly all inflationary models.! Unlike other features of inflation, such as flatness 
or a nearly Harrison-Zel'dovich spectra, the predicted relation between the spectral indices 
and amplitudes was unanticipated prior to the development of the inflationary model and 
is, in this sense, a unique stamp of inflation. 

It is hoped that forthcoming measurements of CMB anisotropy,!' complemented by 
measurements of CMB polarization,^ can be used to detect a tensor component and to 
test the inflationary scenario. The CMB temperature anisotropy power spectrum can be 
expressed in multipole moments where each Cg can be divided into a sum of independent 

rrp\ { a\ 

tensor (C\ ) and scalar (C\ ) subcomponents. While the individual tensor and scalar 
contribution to the temperature anisotropy cannot be observed independently, the presence 
of a strong tensor component may be measured through correlations of polarization and 
temperature.0'E!'0 The amplitude of these correlations, characterized by the ratio of tensor- 
to-scalar multipole moments, = C g jC\ , can be used to test the inflationary scenario. 
However, the moments do not depend on the primordial relation between spectral indices 
and amplitudes alone. The power spectrum also depends on the evolutionary history since 
inflation ended. For example, in the standard cold dark matter (sCDM) model in which 
Vt m = 1, the inflationary predictionHll! is well-known to be 

r 2 « 7(1 -n s ). 

This prediction is modified for open CDM models or CDM models with a cosmological con- 
stant (A), in which the post-inflation evolution is different, as has been studied previously!! 

The purpose of this paper is to discuss the gravitational wave contribution to the CMB 
anisotropy in a new class of cosmological models where a significant fraction of the energy 
density of the universe takes the form of a cosmic scalar field (Q) with an equation-of- 
state different from that of matter, radiation or cosmological constant .0 The scalar field 
component of the cosmic fluid has been dubbed "quintessence" and cosmological models 
based on a combination of quintessence and cold dark matter components are known as 
QCDM models. In a recent paper0 (henceforth, referred to as Paper I), we computed the 
background evolution, the CMB power spectrum and the mass spectrum for QCDM models. 
We showed that models of this type result in a significantly better fit than sCDM to the 
CMB temperature power spectrum, the mass power spectrum, early structure formation, 
and distant supernovae and gravitational lens count measurements. In Paper I we restricted 
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our attention to the case of strictly scale-invariant (ns = 1) spectra of purely scalar metric 
perturbations. In the present paper, we discuss how the tensor contribution to the CMB 
anisotropy is modified by the presence of a Q-component. We compute the full, scalar plus 
tensor CMB power spectrum for QCDM models based on inflationary initial conditions and 
determine how the relation between primordial spectral indices and amplitudes predicted 
by inflation is modified in QCDM models. The result is a generalization of a key test of 
inflationary cosmology. 

The organization of this paper is as follows. In Section |J we briefly review the predictions 
of inflation for the spectrum of initial tensor and scalar perturbations and the resultant 
contributions to the CMB anisotropy. In Section JTJ we discuss how the computation of the 
tensor contribution is modified in QCDM models. (The computation of the scalar spectrum 
has been described in Paper I.) We then present the computed tensor contribution to the 
CMB power spectrum for QCDM models and compare with an sCDM model. In Section 
|TV| we present our key result, the generalization of the inflationary relation between spectral 
index and amplitude for QCDM models. 



II. REVIEW OF INFLATIONARY PREDICTIONS FOR THE SCDM MODEL 

Inflation predicts a nearly scale-invariant power spectrum of energy density (scalar met- 
ric) and gravitational wave (tensor metric) fluctuations. For both subcomponents, the power 
spectrum as a function of Fourier mode k can be expressed to lowest order in terms of the 
Hubble constant and its time-derivatives evaluated when the Fourier mode k was stretched 
beyond the horizon during inflation, at k = aH. Here a is the Friedmann-Robertson- Walker 
scale factor and H = H(<ft) is the Hubble parameter, which depends on the expectation value 
of the inflaton field, 0. Expanding around some given wavenumber ko, the power spectrum 
can be parameterized by spectral amplitudes A(s,t) and spectral indices Ji(s,r) : 

ns-l 

P s (k) = AU-^\ =4 1— — -J \ k=aH (2.1) 
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P T (k) = A 2 T \^-j = — ( — ) \k=aH, (2.2) 

where H' = dH/ dip and m p is the Planck mass. 

A signature of inflation is the series of relations between spectral indices and amplitudes.i 
The tensor spectral index, H(<p) and the power spectrum amplitudes obey the relation§ll! 

1 A 2 

- ~\% < 2 ' 3) 

The ratio of power spectrum amplitudes can be related to the ratio of tensor-to-scalar 
contributions to the CMB anisotropy quadrupole moment 




C (T) 

r 2 = » -7n T . (2.4) 

U 2 
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We will comment on the coefficient below. The scalar spectral index satisfies a more com- 
plicated relation: 




ns _ i = nT _ _^ . (2.5) 



However, for all but an exceptional set of inflaton potentials, the Hubble parameter evolves 
so slowly during inflation that the second term is negligible and ns — 1 ~ nr- In this case, 

r 2 » 7(1 -n s ). (2.6) 

Hence, the relations arise because all parameters associated with the perturbation spectrum 
are ultimately determined only by H and its derivatives during inflation. 

The coefficient in Eqs. ( |2.4| ) and ( |2.6| ) is crucial for the purposes of this paper. In the 
sCDM model, the coefficient is set largely by the Sachs- Wolfe effect, the variation of the 
gravitational potential on the surface of last scattering, which is only very weakly dependent 
on model parameters such as the Hubble constant or the baryon density. Hence, Eq. ( |2.6| ) 
has been presented in the literature!' as a robust prediction of inflation. However, this 
prediction is valid only if Q m = 1. In models in which the matter density is less than unity, 
such as the QCDM models considered in this paper or models with a cosmological constant, 
there is a large, integrated Sachs- Wolfe contribution to the large angle CMB anisotropy, 
which changes the predicted relationship between the spectral indices and the tensor-to- 
scalar ratio. Because the integrated Sachs- Wolfe effect depends on the time variation of the 
gravitational potential along the line-of-site to the last scattering surface, the amplitude of 
the effect is sensitive to the present value of Q m = 1 — Qq, the equation-of-state w, and 
the variation of w with time. Consequently, the coefficient in Eq. ( |2.6| ) is modified by a 
model-dependent function (Section |V|), the central result of this paper. 



III. TENSOR CONTRIBUTION TO THE CMB ANISOTROPY IN QCDM 

MODELS 

In this section we discuss the QCDM models for which we have computed the gravita- 
tional wave contribution to the CMB anisotropy. We present the equations necessary to 
evolve the background and tensor perturbations. We discuss the properties of the Q-matter 
and their dependence on the effective potential for Q, V(Q), and identify two broad cat- 
egories of models. This classification simplifies our survey, in the following section, of the 
imprint of gravitational waves in QCDM models. Finally, we present some sample results of 
our computation of the tensor contribution to the CMB anisotropy spectrum. 



A. Background Equations 

The QCDM models are constructed from spatially flat, Friedmann-Robertson- Walker 
(FRW) space-times containing baryons, cold dark matter, neutrinos, radiation and a cosmic 
scalar field or Q-component. The space-time metric is given by ds 2 = —dt 2 + a 2 (t)dx 2 where 
a is the expansion scale factor and t is the cosmological time. The background equation of 
motion, the energy density, and the pressure for the Q-matter are 
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dV 

Q + 3HQ = - — , (3.1) 
PQ = \Q 2 +V, Pq = \Q 2 - V, (3.2) 



where the dot represents d/dt. Hence, the relations in ( p.l[]3.2| ) supplement the usual back- 
ground equations to specify the evolution of all components of the cosmological fluid. 



B. Gravitational Wave Equations 

Fluctuations about the background space-time can be represented in the synchronous 
gauge, where the linearized space-time metric is 

ds 2 = -dt 2 + a 2 (t)(7i j + h^drfdx?. (3.3) 

Here, 7^ is the unperturbed spatial metric and hij = hfj + hfj is the metric perturbation, 
which includes scalar (S) and tensor (T) perturbations. The equations of motion for the 
scalar perturbations and the fluctuations in the Q-component have been described elsewhere. 
The transverse, traceless gauge constraints and equations of motion of the gravitational wave 
perturbation are 

W« = i ij h T iA = 0, 



d 2 d 1 

— + 3H— V m V 

dt 2 dt a 2 



in 



h 1 l j = 0, (3.4) 



where Vj denotes the covariant derivative with respect to the spatial metric 7^. Because we 
have modeled the Q-matter as a scalar field, the explicit form of the equation of motion for 
the tensor perturbations is unchanged; there are no additional inhomogeneous tensor fields 
or anisotropic tensor sources (as would occur if the Q-component were modeled as a tangled 
web of non-intercommuting cosmic strings^ for example). Hence, the only effect on the 
evolution of the gravitational wave amplitude is through the background expansion. It is 
then straightforward to adopt the standard algorithms for evolving the Boltzmann equations 
with tensor fluctuations for the problem at hand. We have modified a series of standard 
Boltzmann CMB codes based on the synchronous gauge, maintaining a fixed relationship 
between the initial amplitudes of the scalar and tensor power spectra, for the computations 



described below J 
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C. Classification of QCDM Models 

The CMB anisotropy spectrum, both scalar and tensor subcomponents, is sensitive to 
the detailed time evolution of the Hubble parameter H since last scattering. In the case of 
the QCDM models, these subcomponents are affected by the time evolution of the equation- 
of-state of the Q-component, w = Pq/pq, which in turn reflects the form of the potential, 
V(Q). We find that for a wide range of potentials, however, there are two broad categories 
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by which we may classify the behavior of the models. Hence, we will focus our attention on 
representative models from each of these categories. 

The two categories correspond to cases where w increases monotonically versus cases 
where w has begun to oscillate about an asymptotic value by the present epoch. In either 
case, the initial value of the scalar field Q is assumed to be set at a particular value by some 
conditions (e.g., inflation) in the early universe. So long as V"(Q) <C H 2 , Q remains fixed 
because the Hubble red shift term (3HQ) dominates the equation-of-motion. Consequently, 
the energy density is nearly constant during the early history of the universe so that w ~ — 1. 
Once H reduces sufficiently that Q begins to evolve down the potential, the balance of kinetic 
and potential energy in Q changes and w begins to grow. Depending on V(Q) and the initial 
conditions, w may continue to grow up to the present epoch or it may begin to oscillate 
around some asymptotic value. The first possibility corresponds to the "monotonic" class 
of models. Here we find that the imprint of Q on the CMB anisotropy and the mass power 
spectrum is well-approximated by the result obtained if w is held constant at roughly the 
mean value during the period when Qq is non-negligible. See Figure [l|. The oscillatory 
class of models corresponds to cases where Q evolves to a point in its potential where w 
begins to oscillate around some asymptotic value. One example is a potential in which Q 
evolves towards, and then begins to oscillate about, a minimum of V(Q) before the present 
epoch. Correspondingly, w begins to oscillate around a mean value determined by the shape 
of the potential, e.g. w — ► for a quadratic potential (V oc Q 2 ), w — > 1/3 for a quartic 
potential (V oc Q A ). See Figure [| Another example is an exponential potential of the form 
V = m 4 exp(—/3Q); as the potential energy begins to dominate the universe, the equation- 
of-state first overshoots, then relaxes towards w — > (j3 2 /247r) — 1 for (3 < \/487r. We will use 
this exponential example to represent the oscillatory-?/; category. 

Note that a given potential can belong to either the monotonic or oscillatory class de- 
pending on parameters and the initial conditions for Q. Namely, for some choices of the 
initial value of Q, w may have begun to oscillate around its asymptotic value by the present 
epoch. For other choices, Q may not have evolved so far by the present epoch and w has 
been increasing monotonically. For example, while an exponential potential is illustrated 
both in Fig. p] (under monotonic) and Fig. |2] (under oscillatory), different initial conditions, 
and values of m and (3, have been selected in the two cases. 

D. Scalar and Tensor Anisotropy Spectra in QCDM Models 

We have computed the CMB anisotropy spectrum due to the scalar and tensor subcom- 
ponents for a number of QCDM models, which we now present. 

Figure ^| shows the full CMB power spectrum predicted by inflation for tilted sCDM 
and a QCDM model with a constant w = —1/3 equation-of-state and Qq = 0.7, each 
with n$ = 0.9. In both panels the total spectrum is broken down into scalar and tensor 
components. It it noticeable that in the range 2 < i < 10 the fractional contribution of the 
tensor spectrum to the total power is different in the two models. 

Figure |] indicates how the Q-component changes the shape of the tensor subcomponent 
of the power spectrum depending on w for a set of constant equation-of-state models. The 
tensor power spectra in the first panel have been artificially normalized to C\ — 1 to 
compare the shapes, demonstrating that the shape is not strongly affected by the change in 
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w. In the second panel, the curves have been properly normalized with respect to COBE. 
Here we see that the main difference is in the overall amplitude; for increasing w, the 
fractional contribution of the tensor spectrum decreases. In Fig. |5|, a similar set of panels 
demonstrates the effect of changes in Qq on the tensor subcomponent. While the shape is 
not strongly affected, we see that for increasing Qq, the fractional contribution of the tensor 
spectrum decreases. 

In Paper I, it was already noted that, even for purely scalar metric fluctuations, the CMB 
power spectrum at large angular scales (low £) in QCDM models exhibits unusual features 
that do not occur in sCDM or other conventional models. This is owing to a combination of 
the modification of cosmic expansion caused by Q (that is, an integrated Sachs- Wolfe effect) 
and the direct effect of fluctuations in Q. Adding a tensor component can add to further 
features at low I. 



IV. GENERALIZATION OF THE INFLATIONARY PREDICTION FOR QCDM 

MODELS 

We have discussed in Section |TTT| how a very large class of QCDM models can be divided 
into two categories: (1) models in which w is constant or monotonically increasing; and, 
(2) models in which w overshoots and then approaches or oscillates about an asymptotic 
value. We have computed the scalar and tensor components of the CMB anisotropy for 
representative models of each type for a wide range of parameters. We have used the 
numerical results to obtain a revision of the inflationary relation between spectral amplitudes 
and spectral indices. 

Our results are expressed in terms of an empirical relation between the scalar spectral 
index n s and the ratios r 2 = C^/df^ and r 10 = C^/C^p . The quadrupole is a con- 
ventional choice; we have also chosen I = 10 because C\$ is only weakly-dependent on the 
cosmological model (compared to C 2 ) when the predicted spectra are COBE-normalized, 
and so our relations can be applied more simply. These differences between QCDM and 
sCDM results have been expressed in terms of correction factors, f 2 and /10, multiplying the 
known sCDM relationship 



r 2 
rio 



7(1 -n s )] x f 2 (Q Q ,n s ) 

4.8(1 -ns)] xf w (Q Q ,n s ). (4.1) 



We have not included the higher order corrections, such as those proportional to dns/dlnk, 
which are negligibly small for most models.il The functions ^2,10 ar e defined in the following 
subsections. In all cases, / — > 1 as Qq — > 0. The dependence on cosmological parameters 
h, Qb is very weak (which is why the relations are considered to be model-independent tests 
of inflation) and has been ignored. 



A. Monotonic Evolution of w 



Models in which w evolves monotonically leave an imprint on the CMB that is well- 
approximated by the constant w models in which w is set to the average value during which 
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fig is non-negligible. The ratios r 2 and rio for the constant equation-of-state models are 
shown in Fig. ^for the case ng = 0.9. Based on plots of this sort for a range of n$, we have 
obtained the following correction factors for QCDM models. For strictly constant w models, 
empirical fits are valid to within 10% for < w < —1, 0.7 < Jig < 1, and < Qq < 0.7. 



l + i2(l-n s )(2 + u;)ft 



g 2 (n Q /x) 



<? 2 (y) = -0.21 + 2.35y-1.03y 2 , 



x 



/ 10 = 1+ (l- ns )((8 + 7w)Q 2 Q + 3) (l-Q Q /x) 



10 

Pio(y) =0.18 + 0.8% 2 



X 



-w H — w z 
3 3 



(4.2) 



(4.3) 



This somewhat complicated expression is needed to have a formula that fits the dependence 
of all three parameters (w, Qq, and ns) over the stated range. In the case of ACDM, where 
w = — 1, the function f 2 agrees with the ACDM result .111 

B. Oscillatory Evolution of w 



Models in which w overshoots and then approaches or oscillates about an asymptotic 
value are well-represented by exponential potentials in which V(Q) dominates the energy 
density of the universe by the present epoch. The ratios r 2 and r 10 for the exponential 
potential models are shown in Fig. ^. Empirical fits to the correction factors, for the same 
range of models and with the same accuracy as described in the previous subsection, are 
given below. 



h 



1 + H ( i _ ns)i 2 + w) tf Q ] (1 - n Q /xr^ 



92{y) 



-0.21 + 2.35y- 1.03^ 



x — 1 w 

5 
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1 



1 + -(1 - n 5 )((8 + 7w)n 2 Q + 3)J (1 - nQ/x)* * 1 *™ 



g w (y) = 0.18 + 0.84y 2 



x = 1 — -w 



(4.4) 



(4.5) 



In these expressions, w represents the present value of the equation-of-state, w = w(t ). 
In all cases, the initial conditions correspond to w — > — 1. In the expressions above, the 
limit w = — 1 corresponds to w = — 1 throughout, which coincides with the standard ACDM 
result .0 We note that these equations do not apply to extreme cases in which the oscillations 
in the Q-matter are strong enough and begin recently enough to leave a distinct feature at 
large angular scales in the CMB power spectrum (such as a sharp peak at low £). Because 
the scalar and tensor subcomponents for such a model depend sensitively on the detailed 
evolution of the Q-component, a correction factor must be computed model-by-model. 
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V. CONCLUSION 



We have described how Boltzmann codes to compute the CMB anisotropy power spec- 
trum and the mass power spectrum in QCDM models can be simply modified to incorporate 
the contribution of tensor metric fluctuations, as predicted by inflationary cosmology. We 
have demonstrated that a Q-component has two important effects on the tensor component. 
First, by modifying the expansion history of the universe and, hence, producing an inte- 
grated Sachs- Wolfe contribution, a Q-component changes the shape of the tensor anisotropy 
power spectrum (see Fig. |3|). Secondly, since the same effect modifies the scalar component, 
but by a different factor, the ratio of tensor-to-scalar contributions to the CMB anisotropy 
is changed. The net result for any given tilt is to reduce the tensor contribution compared 
to sCDM when spectra are COBE normalized (see Fig. ^). Since the mass power spectrum 
is normalized by the scalar contribution to the CMB anisotropy, a consequence is that, for 
a given tilt, including the tensor contribution does not reduce the COBE normalization of 
the mass power spectrum as much as in sCDM models. (However, as pointed out in Paper 
I, the shape of the mass power spectrum in QCDM is changed and the small-scale power is 
reduced by other effects that do not occur in sCDM models.) 

Finally, the inflationary relations linking the ratio of tensor-to-scalar multipole moments 
to the tilt (jis) must be modified. The most important results presented here are the 
generalized relations shown in the previous section. These relations provide the key test for 
inflation in future CMB anisotropy and polarization measurements, now extended to include 
QCDM models. Assuming inflation is correct, the relations are also used in the fitting 
procedure to determine the Hubble parameter, baryon density, and other cosmic parameters 
from the CMB anisotropyS The determination of cosmic parameters from CMB anisotropy 
measurements in QCDM models, including the effect of the new relations between spectral 
amplitudes and tilt, will be discussed in a forthcoming publication. 
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FIGURES 




FIG. 1. For the monotonic class of QCDM models discussed in the text, the evolution of the 
equation-of-state, w, and energy density, Qq, are shown in the left and right panels. The examples 
shown are for a constant w = —2/3 model and exponential, quartic, and quadratic potentials which 
reach w(to) = —1/3 and = 0.6 by the present day. 



11 




FIG. 2. For the oscillatory class of QCDM models discussed in the text, the evolution of the 
equation-of-state, w, and energy density, Qq, are shown in the left and right panels, for oscillatory 
quadratic and overshooting exponential potentials. In each case, CIq « 0.2 at the present day. 



12 




multipole moment: / multipole moment: / 

FIG. 3. The CMB anisotropy spectrum, decomposed into the scalar and tensor subcomponents, 
for a tilted sCDM model and a tilted QCDM model with a constant equation-of-state uu = —1/3 
and Qq = 0.7. In each case, ns = 0.9, and the relationship between the primordial power spectra 
as predicted by inflation has been maintained. 
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FIG. 4. The tensor subcomponent of the CMB anisotropy spectrum for a series of three QCDM 
models with fig = 0.5 and ns = 0.9. In the left panel, the tensor quadrupoles have been artifi- 

(T) 

cially set to C\ — 1 in order to compare the shapes. In the right panel, the COBE- normalized 
amplitudes have been restored. 




FIG. 5. The tensor subcomponent of the CMB anisotropy spectrum for a pair of QCDM models 
with uu = —1/2, and tilted sCDM, all with ns = 0.9. In the left panel, the tensor quadrupoles 

(T) 

have been artificially set to = 1 in order to compare the shapes. In the right panel, the 

COBE-normalized amplitudes have been restored. 
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energy density: £l Q energy density: £2 Q 



FIG. 6. The ratios T2 and rio for a series of QCDM models with constant equation-of-state w 
and spectral index n$ = 0.9, as a function of VLq. 
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FIG. 7. The ratios r2 and no for a series of monotonic exponential potential QCDM models 
with spectral index ns = 0.9, as a function of SIq. For these models, the equation-of-state w is 
time- varying and w(to) is the present value of the equation-of-state. 
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